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One of the most challenging problems in understanding the structural phase transformations 
in Pu is to determine the energetically favored, continuous atomic pathways from one crystal 
symmetry to another. This problem involves enumerating candidate pathways and studying 
their energetics to garner insight into instabilities and energy barriers. The purpose of this 
work is to investigate the energetics of two transformation pathways for the 5 — > a' trans- 
formation in Pu that were recently proposed pj on the basis of symmetry. These pathways 
require the presence of either an intermediate hexagonal closed-packed (hep) structure or a 
simple hexagonal (sh) structure. A subgroup of the parent fee and the intermediate hexagonal 
structure, which has trigonal symmetry, facilitates the transformation to the intermediate hep 
or sh structure. Phonons then break the translational symmetry from the intermediate hep 
or sh structure to the final monoclinic symmetry of the a' structure. We perform simulations 
using the modified embedded atom method (MEAM) for Pu to investigate these candidate 
pathways. Our main conclusion is that the path via hep is energetically favored and the vol- 
ume change for both pathways essentially occurs in the second step of the transformation, 
i.e. from the intermediate sh or hep to the monoclinic structure. Our work also highlights 
the deficiency of the current state-of-the-art MEAM potential in capturing the anisotropy 
associated with the lower symmetry monoclinic structure. 

Keywords: plutonium, phase transformation, intermediate phase, modified embedded atom 
method 



1. Introduction 

One of the outstanding problems in Pu science is understanding the mechanism of 
the structural transformation from the high-temperature fee 5 phase to the low- 
temperature monoclinic a phase [l|, Q, 0, This transformation is accompanied 
by changes in the shape of the unit cell and numbers of atoms within the unit cells. 
Thus, a combination of strains and shuffles (displacement modes) are involved 
in describing these transformations and the challenging problem is to obtain the 
atomic pathways from one crystal symmetry to another. This is a highly nontrivial 
problem as the nature of the transformation mechanism itself depends crucially 
on orientation relationships between the structures and, therefore, its elucidation 
requires exploring a myriad of symmetry relationships between them. For exam- 
ple, a reconstructive transformation typically involves the presence of intermediate 
structures not related to the parent and the product phases by group-subgroup 
relations that allow the atoms to move continuously, subjected to the constraints 
imposed by the orientation relationships. In contrast, a displacive transformation is 
a one-step group-subgroup process occurring without the need to invoke intermedi- 
ate structures and involves relatively straightforward orientation relations. Recent 
advances in exploring and searching through large databases of crystal symmetries 
0] have made this problem more tractable and provide the means of enumer- 
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ating candidate pathways the energetics of which can be subsequently studied by 
electronic structure or molecular dynamics calculations. 

Lookman et al. [H have recently suggested that the transformation from the fee 
5 phase stabilized by certain trivalent impurities, such as Ga, Al, Ce or Am, to 
the monoclinic a' phase can be accomplished using specific strains, displacement 
modes, and via intermediate structures. This model is consistent with experimental 
data, notably the results of phonon dispersion experiments obtained by inelastic X- 
ray scattering on a thin polycrystalline sample of Pu-Ga alloy [t^. The presence of 
these impurities merely changes the lattice parameters and local ordering but does 
not affect the global crystal symmetry, hence 6 and (Ga, Al, Ce, Am)-stabilized 5, 
as well as a and a', are assumed to have the same respective crystal symmetries. 
On the basis of symmetry, the experimental orientation relationships, i.e. [lll]^ || 
[010]„ and [lIO]^ || [112] Q., impose severe constraints on the 6 — > a{oi') transforma- 
tion. As a consequence, the required intermediate structure turns out to possess 
hexagonal symmetry and the transformation proceeds in two steps. The first step 
is from the fee structure to the hexagonal symmetry facilitated by their common 
subgroup, the trigonal symmetry. The second is from the hexagonal to a monoclinic 
structure, which is a group-subgroup transformation driven by collective displace- 
ment modes, or phonons, that break translational symmetry by combining either 
8 hexagonal close-packed (hep) or 16 simple hexagonal (sh) primitive cells. Fig. [1] 
shows both transformation pathways via the common trigonal subgroup. The con- 
clusion that the intermediate crystal structure has to have hexagonal symmetry 
can be appreciated by recognizing that the only way an fee threefold [111] axis 
can become a monoclinic twofold [010] axis is by an increase of symmetry via a 
structure with a sixfold axis (product of the individual symmetries of the fee [111] 
and the monoclinic [010] axes), that is characteristic for the hexagonal symmetry. 
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Figure 1. The phonon mechanism for (5 — > a (a') phase transformation as proposed in Ref. [ij. The arrows 
indicate displacive (martensitic) phase transformations between structures that are group-subgroup related. 
The labels attached to individual arrows indicate specific phonons and elastic constants associated with 
the strains of the parent structure that drive the corresponding transformation. The Wyckoff symbols of 
each structure are attached to the left of each box. 



Although symmetry is immensely useful in constraining the enormous phase 
space of possible pathways, it needs to be complemented by energy calculations to 
indicate the presence of instabilities and barriers in order to determine the likely 
atomic pathways. First-principles electronic structure and atomistic calculations 
of strongly correlated materials, in particular Pu and its alloys, have long been 
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at the forefront of computational approaches to complement scarcely available ex- 
perimental data and to provide insight into the fundamental physics. The density 
functional theory (DFT) in the usual local density approximation (LDA) or the 
generalized gradient approximation (GGA) predicts a 30% larger volume of the 6 
phase and magnetic long-range order. Electronic structure calculations of Pu based 
on the full-potential linear augmented plane-wave (FP-LAPW) method utilizing 
an antiferromagnetic configuration have shown to reproduce the bulk modulus 
of the a, 7 and 5 phases better than a non-spin-polarized calculation. Currently the 
only first principles approach that has been demonstrated to reproduce all solid- 
state phases of Pu utilizes the full-potential linear muffin-tin orbital (FP-LMTO) 
method with spin-orbit coupling and orbital polarization effects [s!]. In this work, 
antiferromagnetic optimized configurations were used for the a, (3 and 7 phases, 
while an approximation using disordered magnetic moments was adopted for 5, 5' 
and e phases. At present, experiments suggest that there are no static or dynamic 
magnetic moments in the a or 5 phases of Pu [!(?], which may be due to complex 
screening effects that obscure the experimental observations [9|]. The most recent 
method of incorporating strong correlations arising from the almost half-filled 5/ 



orbitals is based on augmenting DFT with the Hubbard Coulomb correlations 
between electrons, giving rise to the "DFT-I-U" method. However, even with this 
additional complexity, the many-body wavefunctions for a crystalline solid do not 
reduce to the corresponding wavefunctions for isolated atoms as the lattice spacing 
is increased [13] . This drawback is removed in the more recent dynamical mean field 
theory (DMFT) in which the Anderson irnpurity model [13] is utilized to study the 



tendency of the / electrons to delocalize |12l ]. 



Despite the success of DFT calculations at capturing correctly the internal ener- 
gies of all crystal structures [§] and the ability of DMFT to reproduce the phonon 
properties of the 5 phase of Pu they are so far impractical as a means to deter- 
mine energies for transformation pathways between the 6 and a{a') phases. These 
simulations could be performed efficiently utilizing a well-parametrized empirical 
or semi-empirical interatomic potential that would capture at least qualitatively 
the strong directional bonding arising from itinerant behavior of the 5/ electrons 
in the low-symmetry monoclinic a and (3 phases of Pu, and their localization in the 
structures of high symmetry such as the orthorhombic 7, face-centered cubic (fee) 6 
or body-centered cubic (bcc) e structures. The current state-of-the-art among these 
approaches is the modified embedded atom method (MEAM), developed specifi- 
cally for Pu by Baskes [3], in which each atom is embedded in the background 
electron density formed by the s, p, d, and / electron densities centered on all other 
atoms. 

The purpose of this work is to formulate a simplified approach that allows stud- 
ies of the energetics of the two transformation pathways proposed in [it] for the 
5 — > a', i.e. fee to monoclinic, transformation via the intermediate hep or simple 
hexagonal structures. We utilize atomistic simulations within the MEAM to probe 
the two candidate paths and calculate the energetics associated with these paths. 
Our principal finding is that the 5 — > a{a') pathway via the intermediate hep struc- 
ture has lower energy barrier than that obtained for the intermediate sh structure. 
Moreover, for both cases the volume change principally occurs in the second step, 
that is, during the group-subgroup displacive transformation from the hexagonal 
to the monoclinic symmetry. Our findings also point to serious deficiencies in the 
current state-of-the-art MEAM potential for Pu which are presumably due to its 
inability to reproduce the experimental positions of atoms in the low-symmetry 
a and (3 phases of Pu. We conclude by outlining possible improvements of these 
calculations that may allow for future systematic search for transformation path- 
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ways and energetics without a priori specifying the space group of the intermediate 
crystal structure. 

2. Interatomic potential and equilibrium crystal structures 

In the following calculations, we utilize the modified embedded atom method 
(MEAM) interatomic potential for Pu [15], within which the total energy is written 
as 

Here, Eemb is an energy to embed an isolated atom into a background electron den- 
sity formed by electrons centered on all other atoms, and Epair is a pair potential. 
The embedding energy takes a form that is often used for the exchange-correlation 
functional in the LDA approximation of the DFT. The adjustable parameters in 
the functional form of the embedding energy are calculated so as the potential 
reproduces some fundamental properties of the fee 6 phase, e.g. the experimental 
elastic constants. The pair potential is then constructed so as the dependence of 
the total energy on volume for the reference fee structure matches the universal 
equation of state [l6^ that is parametrized by the experimental bulk modulus of 
the 5 phase and by the pressure derivative of the bulk modulus of the a phas43- If 
the theoretical formalism of this potential captures the essential physical aspects 
of bonding at least qualitatively correctly, one might expect that the properties 
associated with the 5 phase and other high-symmetry structures that do not de- 
viate significantly from this reference lattice will be reproduced as well. On the 
other hand, the monoclinic a phase has 20% lower volume than the reference fee 
structure and the potential may not be capable to reproduce the experimental 
measurements with sufficient accuracy. 

In order for our simulations to be consistent with the interatomic potential used, 
it is necessary to obtain the positions of atoms and the lattice parameters of the 
monoclinic cell that corresponds to the minimum of energy. This is accomplished 
by utilizing the isothermal-isobaric molecular dynamics simulation at K where 
the temperature is maintained constant by simple rescaling of velocities to obey 
the equipartition of kinetic energy and the pressure is controlled by the modularly 
invariant Cleveland- Wentzcovitch barostat [3, [3] ■ The latter allows the simulated 
cell to change shape as the atomic positions are evolved during the minimization of 
energy. The ensuing system of coupled equations of motion for the atomic positions 
and for the shape of the simulated cell (i.e. lattice parameters a, b, c and the lattice 
angles a, f5, 7) is integrated numerically using the fourth order Gear predictor- 
corrector method [13] with constant time step. To speed up the relaxation process 
we constrain the atoms to obey the twofold symmetry of the b axis of the monoclinic 
lattice. In particular, the equations of motion for the atoms are solved only for the 
eight nonequivalent atoms in one of the two monoclinic planes while the positions 
of the other eight atoms in the second plane are calculated from the symmetry 
mapping (a, |,c) <-> (a, |,c), where a and c are fractional coordinates of atoms 
measured relative to the a and c axes. The energy of the relaxed monoclinic lattice 
is about —4.1 eV/atom which is in a good agreement with the values extrapolated 
from experimental data [isl ]. Similarly, the energy difference between the 6 and a 
structures used in our simulations is about 0.3 eV/atom, which agrees with the 



^Similar measurements for the S phase were not available at the time. These experiments were performed 
recently and are published in Ref. [l7j . 
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total internal energy differences estimated in [21|] by accounting for the thermal 
expansion of the lattice within each phase. The calculated volume of the a structure 
is 20.8 A'^/atom which slightly overestimates the measured experimental volumes 
extrapolated to K. 



3. Energetics of transformation pathways 

In the most general situation the prominence of individual phonon modes varies 
along the transformation pathway, and this necessarily gives rise to the motion of 
atoms along curvilinear paths. Unless one performs a more sophisticated simula- 
tions, such as the nudged elastic band (NEB) method [i^], which cannot be done 
at present due to apparent difficulties of the MEAM potential to reproduce the 
positions of atoms in the a structure (more on this later), the information about 
the prominence of individual phonon modes is not available. Hence, one is forced 
to assume that the superposition of a number of phonons responsible for a given 
transformation results in the motion of atoms along the shortest paths connecting 
their initial and final positions, which is a simplification that we also adopt in our 
simulations here. Although these calculations were carried out for pure Pu, our 
conclusions are qualitatively valid also for dilute alloys of Pu with Ga, Al, Ce or 
Am that stabilize the fee b phase down to essentially K. The reason is that low 
concentrations of these impurities do not change significantly the lattice parame- 
ters. In particular, additions of Ga and Al slightly shrink the fee lattice while Ce 
and Am cause its expansion Nevertheless, for low concentrations of these al- 
loying elements the changes in the lattice constant do not exceed 1% of the values 
for pure Pu. 

3.1. fee — 'r sh ^ monoclinic transformation 

The simulated block is constructed by placing the atoms in their equilibrium lat- 
tice positions in the fee structure, as shown in Fig. [2ja). Due to the orientation 
relationship [111] 5 || [0001] sh, the minimum number of atoms in the initial fee struc- 
ture is three. During the course of the transformation, these atoms are considered 
to move along the shortest paths from their initial positions, depicted as gray in 
Fig. EJa), to the three corners of the final hexagonal lattice in the directions shown 
in the same figure by arrows. The extent of this transformation is conveniently 
measured by the transformation coordinate x that represents the displacement of 
each atom from its initial (fee) to the final (sh) lattice. Here, x = corresponds 
to the fee lattice, < x < 1 to the intermediate trigonal lattice, and x = 1 to the 
final sh lattice. Each value of x thus corresponds to different fractional coordinates 
of atoms that are calculated as 

fi(a;) =f° + xt; , (2) 

where the hat {') designates fractional coordinates measured relative to the shape 
of the simulated cell, and tj is the vector that displaces the atom i from its ini- 
tial position in the fee structure to its final position in the sh lattice. For known 
positions of atoms, given by fj(x), we calculate the shape of the unit cell that 
minimizes the total energy of the system. During this optimization the fractional 
coordinates of atoms do not change and the absolute motion of atoms is only due 
to the changing shape of the simulated cell. The relaxed configuration is character- 
ized by the lattice parameters a, b, c, the lattice angles a, /3, 7, the corresponding 
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volume of the cell V and the energy E. The variation of the volume and energy as 
a function of the transformation coordinate x is shown in the left panel of Fig. [3l 
The activation barrier for this transformation is about 0.1 eV/atom. An impor- 
tant aspect to note is that no volume change is observed after the transformation 
is completed. If the sh phase were the intermediate phase, it would thus mean 
that the 20% volume decrease during the 5 a{a') transformation observed in 
experiments has to occur during the sh monoclinic transformation. 



fee ^ sh sh — » monocUnic 




(c) (d) 



Figure 2. Motion of atoms during the fee — > sh (a), sh — > monoelinie (b), fee — > hep (c), and hep — * 
monoeUnie (d) transformations. The atoms in the parent lattiee (outlined) are plotted in gray and their 
corresponding positions after the transformation is complete are in black; the shape of the final cell is not 
shown for clarity. The arrows show the motion of the atoms between the two structures, i.e. the vectors ti 
in ig. 

For the calculation of the sh monoclinic transformation we utilized the relaxed 
sh lattice from above to construct a 2 x 2 x 4 sh supercell. The positions of atoms in 
this initial supercell are shown in Fig. E^b) as gray while their final positions in the 
monoclinic lattice are in black. Although the assignment of atomic correspondences 
between the two crystal structures (i.e. sh and monoclinic) is not unique, one can 
expect that the collective motion of the atoms proceeds in such a way that the 
activation energy for this transformation is minimized. The corresponding atomic 
correspondences are shown in Fig. EJ^b) by arrows that designate the direction of 
motion of each atom and thus the vector tj. From the Shoji-Nishiyama relationship, 
the [0001] axis of this sh cell will coincide at the end of the transformation with 
the [010] axis of the final monoclinic structure. This helps us to identify the corre- 
spondences between the positions of atoms in the initial and final crystal structures 
that are determined from the assumption that each atom in the initial sh supercell 
moves to its corresponding position in the monoclinic cell in the shortest distance 
when measured in fractional coordinates. The transformation coordinate x now 
represents the displacement of atoms from their initial position in the sh supercell, 
i.e. X = corresponds to the initial sh lattice and x = 1 to the final monoclinic 
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lattice. The instantaneous position of each atom is again calculated using ([2]) with 
the displacement vectors tj indicated in Fig. Wijo). 
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Figure 3. Variation of volume (V) and energy (E) during the transformation of Pu from its fee 5 phase 
to the monoclinic a phase. In the left panel of the figure the atoms are moving from their initial positions 
in the fee lattice to their final positions in the sh (hep) lattice, while in the right panel of the figure the 
atoms are moving from their positions in the sh (hep) lattice to their final positions in the monoclinic 
lattice. The transformation coordinate x designates the progress of the transformation (x = corresponds 
to the initial structure and x = 1 to the final structure). In the right panel of the figure, the solid lines 
correspond to the 16-atom conventional monoclinic cell, that is stable within MEAM, with angle 106.3° 
and the dashed lines to the 32-atom monoclinic supercell with angle 120.3°. 



The calculated dependence of the relaxed volume of the cell and the correspond- 
ing energy is shown in the right panel of Fig. [3l As expected above, the volume 
necessarily undergoes about 20% decrease. The energy increases by an additional 
0.2 eV/atom and then decreases to the final value of about -4.1 eV/atom which is 
the energy of the monoclinic lattice that is stable within MEAM. Accounting for 
the energy barriers from the two transformation steps, the total energy barrier for 
the fee — > sh — > monoclinic transformation becomes about 0.3 eV/atom. 



3.2. fee — > hep — > monoelinie transformation 

For this simulation the initial block of atoms contains six atoms in the six consec- 
utive (111) planes of the fee lattice. The positions of atoms in the initial fee lattice 
are shown in Fig. [2]^c) in gray while their corresponding positions in the final hep 
lattice are in black; the arrows indicate the direction of motion of each atom. After 
the atoms are displaced in straight paths towards their corresponding positions in 
the hep structure, where the instantaneous positions of atoms are calculated from 
(l2|), one obtains three hep cells stacked on the top of each other in the [0001] direc- 
tion. The calculated dependence of the volume and energy on the transformation 
coordinate x is shown in the left panel of Fig. [3l As in the case of fee — > sh step 
discussed above, there is no volume change after completing this transformation 
and thus it is inevitable that the 20% volume decrease occurs in the second step, 
i.e. during the hep — > monoclinic transformation. Interestingly, there is a very low 
energy barrier for the fee — > hep transformation which implies comparable stability 
of the fee and hep phases of Pu. 

In order to investigate the hep — > monoclinic transformation, we again construct 
the initial hep supercell that now consists of a 2 x 1 x 4 array of the hep unit cells 
with two atoms per cell. The initial positions of atoms are shown in Fig. EJd) in 
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gray and their final positions in the monochnic lattice are in black. The arrows 
indicate the displacement vectors U that are substituted in ([2]) to calculate ri{x), 
the instantaneous fractional coordinates of atoms corresponding to a given value 
of the transformation coordinate x. The evolution of the volume and energy of the 
simulated cell is shown in the right panel of Fig. [3l One again observes the 20% 
volume decrease and an additional energy barrier of about 0.2 eV/atom. However, 
since there is virtually no energy barrier observed in the first step of this transfor- 
mation, the total energy barrier for the fee — > hep monoclinic transformation is 
only about 0.2 eV/atom and thus lower than for the case when the intermediate 
phase has sh structure. 



4. Discussion 



The results of the simulations presented in Fig. [3] imply that the preferred inter- 
mediate phase for the fee (6) — >■ monoclinic (a, a') transformation in Pu is the hep 
phase and the corresponding energy barrier is about 0.2 eV/atom. It is well-known 
that Pu has a unique position among the actinides that stems from the behavior 
of its / electrons |23l. l24l|. In particular, the / electrons in early actinides (Ac to 
Np) are itinerant and participate in bonding, whereas in the late actinides (Pu to 
No) they are localized at individual atoms. Since the late actinides (in particular 
Am to Cf) all crystallize in hep (or double-hcp) structures and owing to the fact 
that Pu technically belongs to the group of late actinides, it is not unreasonable 
to imagine that the hep structure should be an intermediate metastable phase me- 
diating the transformation from the fee to the monoclinic phase. An hep phase 
has been recently observed during the pressure-induced (5 — > a' transformation in 
Pu-2at.% Ga alloys IT]]. Utilizing Rietveld refinement, the crystal structure of this 
intermediate phase was found to belong to the orthorhombic Fddd group with the 
atomic positions corresponding to those in slightly distorted hep planes. This inter- 
mediate phase is presumably short-lived (the duration of its stability is of the order 
of picoseconds) and, therefore, further transient X-ray diffraction experiments are 
needed to confirm its existence. 

When investigating the hep monoclinic transformation path, we assumed that 
a 16-atom sh (hep) cell transforms to the conventional monoclinic unit cell with 
16 atoms and the monoclinic angle 106.3°. However, the choice of the tiling of 
the monoclinic lattice is not unique and one can construct an infinite number of 
supercells that represent it, although with larger number of atoms. For example, it 
has been shown [25| that one can construct a supercell with 32 atoms (see Fig. S]) 
whose monoclinic angle is, in our case, 120.3°. This angle is very close to that of 
the initial hexagonal (sh or hep) lattice and, therefore, basically no strain is as- 
sociated with changing the shape of the simulated cell during the transformation. 
One would thus expect that using this 32-atom monoclinic supercell, the energy 
barrier for the sh (hep) — > monoclinic transformation would be lower than that for 
the conventional 16-atom monoclinic cell. We performed the simulations above also 
using this larger 32-atom cell and observed only minor changes in the variation of 
energy and volume for the sh — > monoclinic transformation. This is shown by com- 
paring the solid and dashed lines in the right panel of Fig. [3] that correspond to the 
16-atom and 32-atom cells, respectively. In both cases of the 32-atom calculation 
the activation energy decreased relative to the equivalent 16-atom simulation by 
an additional 0.05 eV/atom. This brings the total energy barrier for the fee — > sh 
— > monoclinic transformation down to about 0.25 eV/atom and that of the fee — > 
hep — > monoclinic transformation to 0.15 eV/atom. 

In order to understand the limitations of the calculations presented in this paper. 
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Figure 4. Four conventional monoclinic cells with 16 atoms and an alternative supercell with 32 atoms 
(filled) both of which represent the monoclinic crystal. The angles of the two cells are 106.3° (16-atom 
cell) and 120.3° (32-atom cell). The atoms in the two monoclinic planes are drawn as black and gray; the 
empty circles correspond to the atoms that are present in the four 16-atom cells but lie outside of the 
32-atom cell. 



it is important to analyze the accuracy with which the MEAM potential represents 
the individual solid phases of Pu. The reference phase that is used to construct 
this potential is the fee 5 phase for which the elastic constants [2^ . 27, H^, lattice 



parameters and the cohesive energy jl5l | are well known. The theoretical formalism 
of any interatomic potential captures only certain relevant aspects of the physics of 
bonding. The potential is typically constructed to reproduce experimental or first 
principles data for a simple crystal structure (often cubic) and is tested for transfer- 
ability to other atomic environments. In the case of the MEAM potential for Pu, one 
assumes that parametrizing this potential to the d phase will provide an empirical 
description that will reproduce correctly the phases whose symmetry and volume 
are comparable with those of the reference structure. Indeed, it has been shown 
in [3] that this potential also predicts correct energies of the orthorhombic 7 and 
bcc e phases. It is, however, much more challenging to achieve the same degree of 
accuracy for the low-symmetry simple monoclinic a and body-centered monoclinic 
P phases whose volumes are 20% and 10% lower, respectively, than the reference 
fee 5 phase. To demonstrate this we show in Fig. [5^a) and (b) a comparison of the 
equilibrium monoclinic lattice calculated from MEAM and used throughout this 
paper with that obtained by Rietveld refinement of the X-ray powder diffraction 
measurements j29( |, respectively. Clearly, the positions of atoms calculated from the 
MEAM potental for Pu do not correspond to the experiment. The most notable 
feature is an unphysical hexagonal symmetry that can be seen clearly by rotating 
the unit cell about the horizontal a axis. The origin of this discrepancy is presum- 
ably the inability of the potential to account for strong directional bonding arising 
from the overlap of the / orbitals as the volume of the crystal is decreased and 
the symmetry of the parent fee phase (48 symmetry operations) is broken into 12 
variants of the monoclinic lattice (4 symmetry operations). 

The observation that the current MEAM potential for Pu does not represent 
with reasonable accuracy the positions of atoms in the a phase of Pu raises serious 
concerns when calculating the transformation pathway between the 6 and a phases. 
The region where MEAM is least accurate, and thus the results of our simplified 
studies presented above are least trustworthy, corresponds to the shaded area in 
Fig. [3l For this range of transformation coordinates x, the volume and thus also the 
lattice parameters (not shown here) change rapidly as the sh (hep) phase transforms 
to the trigonal phase and then to the final monoclinic structure. Nevertheless, the 
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a a 



(a) (b) 

Figure 5. Atomic positions in the monoclinic a structure calculated using MEAM (a) and those deduced 
from the X-ray powder diffraction measurements of 1291 (b) . The same colors of atoms belong to the same 
monoclinic plane in the b direction that is perpendicular to the plane of the figure. If the lattice in (a) is 
rotated around the horizontal a axis, the black and gray atoms in the two monoclinic planes overlap and 
one observes a perfect hexagonal pattern that is not present in (b). 

inaccuracy of the MEAM potential in this low-volume region does not affect our 
main conclusion that the hep structure is predicted as the preferred intermediate 
symmetry through which the 5 phase transforms to the monoclinic a phase. 



5. Conclusions 

Our objective in this paper was to investigate two prototypical transformation 
mechanisms for Pu that were identified in [l| from purely symmetry considera- 
tions. Utilizing the well-known Shoji-Nishiyama relationship [lll]^ || [0001]hex || 
[010]a that constrains the mutual orientations between the individual phases of 
Pu, Lookman et al. [l[ suggested that the transformation of the fee 5 phase to 
the monoclinic a phase proceeds via the intermediate trigonal and either sh or 
hep structures. However, since these considerations do not involve any energetics, 
they cannot determine which of these hexagonal structures is more energetically 
favorable. 

In the present work, these calculations were performed using the MEAM for Pu 



15| in which we assumed that the fractional coordinates of atoms change along the 
shortest paths between the initial and the final crystal structures. The individual 
phonon modes along the transformation pathway correspond to displacements of 
atoms in different directions but the vector sum of these individual vectors yields 
a straight path between the positions of atoms in their parent and product phases. 
The corresponding shape of the unit cell that minimizes the total energy is deter- 
mined while keeping the fractional coordinates of atoms fixed. These calculations 
are performed in the isothermal-isobaric ensemble at K. We have shown that the 
current MEAM potential for Pu predicts that the 5 a{a') transformation occurs 
preferably via the intermediate hep phase. If one uses the 32-atom supercell for the 
hep — > monoclinic transformation, the total energy barrier is about 0.15 eV/atom 
and the 20% volume decrease occurs entirely during the transformation from the 
intermediate hep to the final monoclinic structure. The existence of an interme- 
diate hep phase has been observed recently during the pressure-induced 6 ^ a' 
phase transformation For completeness, the signatures of the intermediate 
trigonal or hep symmetries should also be probed with picosecond resolution by, 
for example, transient X-ray diffraction techniques that are becoming increasingly 
important in many areas of materials science. 
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We argue that while the MEAM potential for Pu reproduces faithfully the crystal 
structures whose symmetries and volumes are comparable with those of the refer- 
ence fee 6 phase, its current formulation cannot capture the positions of atoms in 
the low-symmetry a and /? phases of Pu that were determined earlier by Rietveld 
refinement of X-ray powder diffraction data i30s]. Hence, also the energetics of 
the pathways from the hexagonal intermediate symmetry to the final monoclinic 
structure is not faithfully reproduced by MEAM. This makes even this state-of- 
the-art interatomic potential inappropriate for use in the methods that search for 
minimum energy paths between two well-defined configurations, such as the NEB 
method 

Continued progress in this field relies on our future ability to understand and 
at least qualitatively correctly capture the core physics associated with strong 
correlations among / electrons. Having an interatomic potential that would suc- 
cessfully reproduce all six crystal phases of Pu and the seventh under pressure 
would undoubtably provide new understanding of the physics of this important 
metal. This potential will have to reproduce the largest set of requirements ever 
imposed on an interatomic potential of a single element. In the case of Pu, these 
are mainly the energies and volumes of all crystal structures, softening of the (111) 
transverse acoustic phonon at the L point in the 6 phase that drives the phase 
transformation to the a' phase, and positions of atoms in the low-symmetry a and 
P structures. One promising candidate is the modified generalized pseudopoten- 



tial method (MGPT) developed by Moriarty et al. [31[ and successfully applied 
to Mo and Ta. Recent advances of this method allowed simulations of / electron 
materials such as 5 Pu in which the strong correlations are modeled by turning off 
the bonding due to the / electrons completely and describing the material as a d 
band metal. Unlike the MEAM potential for Pu, this MGPT potential has been 
demonstrated to capture reasonably well the phonon dispersion curves, although 
the anomalous softening of the transverse (111) phonon at the L point is not re- 
produced. Another possible candidate for the construction of a better interatomic 
potential is the tight-binding-based description within the bond order potential 
(BOP) [l^] that has been shown to be imminently suitable for transition metals 
crystallizing in the bcc structure with approximately half-filled d bands (V, Nb, 
Mo, Ta, W) and a Fe. Although this BOP is not currently available, the first step 
was undertaken by Hachiya [33] who constructed a tight-binding description of fee 
a phase of Th and 5 Pu, both of which are characterized by approximately half- 
filled / orbitals. However, more extensive tests of this approach, in particular its 
ability to describe the space group of the monoclinic a phase of Pu and the large 
volume changes between the 5 and a, (3 phases of Pu is still lacking. 
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